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Lossy chaotic electromagnetic reverberation chambers: 
Universal statistical behavior of the vectorial field 
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The effective Hamiltonian formalism is extended to vectorial electromagnetic waves in order to 
describe statistical properties of the field in reverberation chambers. The latter are commonly 
used in electromagnetic compatibility tests. As a first step, the distribution of wave intensities in 
chaotic systems with varying opening in the weak coupling limit for scalar quantum waves is derived 
by means of random matrix theory. In this limit the only parameters are the modal overlap and 
the number of open channels. Using the extended effective Hamiltonian, we describe the intensity 
statistics of the vectorial electromagnetic eigenmodes of lossy reverberation chambers. Finally, 
the typical quantity of interest in such chambers, namely, the distribution of the electromagnetic 
response, is discussed. By determining the distribution of the phase rigidity, describing the coupling 
to the environment, using random matrix numerical data, we find good agreement between the 
theoretical prediction and numerical calculations of the response. 

PACS numbers: 41.20.-q,42.25.Bs,05.45.Mt 


I. INTRODUCTION 

For more than 40 years, wave chaos has been an at¬ 
tractive field of fundamental research concerning a wide 
variety of physical systems such as quantum physics [1], 
room acoustics [2] or ocean acoustics [3], guided-wave 
optics [4], microwave cavities [5, 6], etc. The success of 
wave chaos is mainly due to its ability to describe such a 
variety of complex systems through a unique formalism 
which permits us to derive a universal statistical behav¬ 
ior. Since the early 1990’s, quasi bi-dimensional chaotic 
electromagnetic (EM) cavities have been one of the most 
used experimental set-ups where these statistical proper¬ 
ties have been verified [7, 8]. In a more applied context, 
EM cavities are nowadays currently used to realize elec¬ 
tromagnetic compatibility (EMC) testing on devices with 
embedded electronic components [9]. In the EMC com¬ 
munity, EM cavities are called EM reverberation cham¬ 
bers (RCs). Thanks to mechanical [9-11] or electronic 
stirring [9, 11-13], enabling ensemble averaging, and to 
the presence of losses leading to modal overlap, the sys¬ 
tems under test are submitted to a supposedly statisti¬ 
cally isotropic, uniform, and depolarized electromagnetic 
field. The statistical description of the EM held inside 
RCs commonly used by the EMC community was origi¬ 
nally proposed by Hill [14] , who made a continuous plane 
wave spectrum hypothesis. This hypothesis assumes that 
the held is statistically equivalent to a random superpo¬ 
sition of traveling plane waves. This hypothesis is gen¬ 
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erally well verihed if the excitation frequency is much 
larger than the so-called lowest usable frequency (LUF), 
the latter being commonly considered to lie between three 
to six times the fundamental frequency [15]. However, in 
a frequency regime close to the LUF, the EM held might 
be neither uniform nor isotropic and this even with stir¬ 
ring [16, 17]. In this regime, conventional RCs display a 
highly non universal behavior depending on their geome¬ 
try [15, 18], the kind of antennas used [19], or the object 
under test [20]. 

Recent investigations of RCs have concentrated on two 
goals: (a) improving the RCs’ behavior for a frequency 
regime close to the LUF and (b) hnding a suitable sta¬ 
tistical model for the EM held in this regime [17, 21-23] 
with the incentive of proposing more accurate quantities 
[23] than those proposed by the lEC standard [9] . In the 
case of (a), the EMC community has mainly focused on 
the optimization of the stirring [18, 24-27] or on the new 
design of the geometry of RCs [28-30]. None of these 
approaches are used to propose quantitative hypotheses 
concerning the statistics of the EM held near the LUF. In 
case of (b), the proposed statistical descriptions converge 
to Hill’s hypothesis in the high frequency limit but are 
hardly justihed on physical grounds. 

Here we suggest a model based on random matrix 
theory (RMT) for open chaotic systems to improve the 
present situation. Indeed, electromagnetic cavities, and 
more particularly EM reverberation chambers are gen¬ 
uine examples of open wave systems due, for instance, to 
Ohmic losses or antenna couplings. The openness of such 
cavities turns the real discrete spectrum of closed sys¬ 
tems, associated to real eigenhelds, into a set of complex 
resonances. The latter manifests themselves as complex 
poles of the scattering matrix which are given by the 
eigenvalues of the effective non-Hermitian Hamiltonian 
HqQ of the open system [32-35] . The anti-Hermitian part 
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FIG. 1. (color online) Lossy parallelepipedic reverberation 
chamber with length /=0.985 m, width w=0.785 m, and height 
/i=0.995m. It is made chaotic through the introduction of a 
movable half-sphere with radius rm=0.15m and two spheri¬ 
cal caps with radii of 55 cm and 45cm[31]. Losses are mainly 
localized on sub-surfaces of the walls with a reduced con¬ 
ductivity (colored patches). The total area of these 

sub-surfaces is denoted There are two configurations, 

namely, cavity 1 where only the light red crossed patches 
(3 X patch 1) are attached and cavity 2 where all six patches 
are attached (3xpatchl and 3 x patch 2). 


of arises from coupling between the internal (bound) 
and continuum states, leading to finite resonance widths. 
The other key feature is that the eigenfunctions of 
are nonorthogonal [33, 35, 36]. For systems invariant un¬ 
der time reversal, like the open EM cavities studied be¬ 
low (see Fig. 1), the nonorthogonality manifests itself in 
complex wave functions, yielding the phase rigidity [37- 
39] and mode complexness [36, 40, 41]. In practice, the 
EM response of such open systems consists of a sum of 
overlapping resonances, associated to complex eigenfields 
(the resonant or decaying states) which can be viewed as 
a sum of standing and traveling waves [42]. 

In this paper we will show how the Ffgff approach [8] 
allows to extend the random matrix theory to fully three- 
dimensional (3D) open chaotic EM reverberation cham¬ 
bers. Hence, universal statistical spectral and spatial be¬ 
haviors can be predicted for these systems. This has been 
recently illustrated in Ref. [36], where the experimental 
width shift distribution was investigated in a chaotic re¬ 
verberation chamber as a test of non-orthogonality. We 
will show here that these crucial properties for electro¬ 
magnetic compatibility tests, especially the field unifor¬ 
mity, can indeed be achieved in chaotic EM cavities [31] 
even at low enough frequency where the modal overlap 
is only weak or moderate. In a conventional RC with a 
regular rectangular shape, these low-lying modes display 
highly non-isotropic patterns resulting in a low-frequency 
response which cannot be expected to be statistically 
isotropic and uniform, even in the presence of stirring. 
To the contrary, these statistical requirements are natu¬ 
rally fulfilled by the vast majority of modes in a chaotic 
cavity without the help of any stirring process [31]. This 
generic statistical behavior of modes in a chaotic cav¬ 
ity will be referred to as ergodicity in the following. In 
particular, we will show how the ergodicity of modes in a 


chaotic RC permits to establish a statistically uniform re¬ 
sponse. This behavior is in complete contradiction with 
what can be obtained in a conventional RC where the 
specific details of the chamber lead to unpredictable sta¬ 
tistical features [43, 44]. 

In the following section, we introduce important nota¬ 
tions and remind the reader of the main universal prop¬ 
erties of open chaotic systems, in particular concerning 
the spatial distribution of complex eigenstates. Most 
are obtained through the effective Hamiltonian formal¬ 
ism and have been successfully verified experimentally 
and numerically in quasi-2D lossy chaotic cavities (see, 
for instance. Ref. [8] and references therein). In the third 
section, we use the Dyadic Green’s function representa¬ 
tion of the response of a fully 3D vectorial EM cavity to 
extend the previous results by properly defining the bi¬ 
orthogonality of complex modes for vectorial fields and 
their phase rigidity. Then, we establish the universal 
statistics of the complex modes and of the full response 
in chaotic RCs. We successfully compare all the pre¬ 
dictions deduced from our extended RMT approach to 
numerical simulations of realistic RCs. We conclude by 
showing how a minimal number of ingredients, namely, 
the number and the coupling strength of fictitious equiv¬ 
alent channels, is sufficient to give a full statistical ac¬ 
count of the intensity distributions that can be measured 
in a chaotic RC, irrespective of the specific details of the 
cavity shape and of the antenna set-up. 


II. CHAOTIC OPEN SYSTEM: UNIVERSAL 
STATISTICS WITH FORMALISM 

The response of an open scalar wave system can be 
formally recovered by means of the Green’s functions for¬ 
malism used in the quantum context for open systems, 
through [45] 

G(e) = (e-iL3fj)-i, (1) 

where TLgff is the (non-Hermitian) effective Hamiltonian 
and is commonly written as [35, 45, 46] 

Heff = i7-^W, (2) 

where H is the Hermitian part corresponding to the 
closed system (neglecting the Hermitian effect due to 
opening) and —^W is the anti-Hermitian part related to 
coupling to the environment. In systems with preserved 
time reversal symmetry (which is assumed for EM cavi¬ 
ties [7, 8]), H is a real symmetric operator. We assume 
W to be real as well since any imaginary part induced 
by the coupling can be included in H[8, 33, 47]. The 
complex eigenvalues of HqQ, = £n — *r^/2, yield the 
center e„ of the nth resonance and its width F^. As 
a consequence of the non-Hermiticity of Ffgff, there are 
two distinct sets of complex left {(L„|} and right {|i?„)} 
eigenvectors, associated to the same set of eigenvalues 
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{en}: 

— (:n\Rn^ • (3) 

The eigenvectors (L„| and \Rn), which describe the res¬ 
onance states, form a complete and bi-orthogonal set: 

Enl^»)(T„| = l, (4) 

{Ln\Rm) — ^nm- (b) 

The time reversal symmetry assumption adds another 

constraint on right and left eigenvectors [33, 35, 41], 
namely that they are the transpose vector of each other: 

(L„| = (6) 

Note that the bi-orthogonality condition (5) implies that 
the real and imaginary parts of the wave functions are 
orthogonal to one another. The Green’s function (1) can 
be expanded on the left and right eigenvectors [35, 46]: 

(7) 

n ^ 

For the present study, important quantities are the mean 
level spacing A and the local average resonance width 
T'^, whose ratio permits to evaluate the amount of over¬ 
lap between adjacent resonances through the mean modal 
overlap: 


A. T^efF model for chaotic cavity and RMT results 

RMT is a powerful tool to derive universal statistical 
properties of closed [7, 48, 49] and open wave chaotic 
systems [8, 40, 50-54]. In this framework, the effective 
Hamiltonian of an open cavity with N resonances is mod¬ 
elled by an TV X A random matrix and loss mechanisms 
are associated with M open channels connecting the N 
levels of the closed (lossless) cavity to its environment. 
As we have in mind to study the statistics of the vec¬ 
torial field of actual chaotic reverberation chambers, we 
concentrate on systems which preserve time reversal sym¬ 
metry. Then the Hamiltonian of the closed system, H in 
Eq. (2), is modelled by a random matrix belonging to 
the Gaussian Orthogonal Ensemble (GOE) [7, 48]. The 
anti-Hermitian part — |IF of is related to a. N x M 
coupling matrix V through W = . Each element 

of V connects the n-th eigenstate of iV {n G [1, A^j) to the 
j-th open channel (j G [1,M]). In order to preserve or¬ 
thogonal invariance of TVgff under orthogonal transforma¬ 
tions [33], the are commonly set as real independent 
Gaussian random variables with zero mean and standard 
deviation aj depending only on channel j (as a conse¬ 
quence of the ergodicity of eigenstates) and related to its 
coupling strength kj through [33, 40, 41] 

cr| = (2KjA)/7r. (9) 


In the following, we will consider equivalent channels as¬ 
sociated to Ohmic losses and will therefore assume that 
all channels are identically coupled [40, 41]: 

(V^Vp^) = ( 10 ) 

We will now establish a few results, obtained in the 
framework of RMT, concerning the statistical features 
of scalar eigenmodes for chaotic open wave systems, as 
verified through experiments in 2D chaotic microwaves 
cavities, and see, in Sec. HI, how they can be extended 
to the case of vectorial fields. 


B. Spatial distribution of the intensity of 
eigenstates 


For scalar fields, eigenmodes of open chaotic systems 
are complex wave functions = (r whose spa¬ 

tial statistics can be described by the statistics of the N 
components (i = 1, • • ■, AT) of the right eigenvectors 
of the random matrix The bi-orthogonality condi¬ 

tion (5) and GOE assumption imply: 


N 


Re [4'„,j] Im [^-n,*] = 0, 

(11) 

N 

(Rc [4'n,i]^ - Im [4'n,z]^) = 1, 

2^1 

(12) 

N 

(Ue + Im [4'n,i]^) > I, 

(13) 


i=l 


where Re [d^n j] and Im [T^ j] are two independent Gaus¬ 
sian distributed random variables. It is important to note 
that Re [dT^^i] and Im [Tn^i] do not necessarily have the 
same variances. The complexness parameter, orig¬ 
inally introduced by Lobkis and Weaver in Ref. [55], is 
defined by the ratio of the variances of the imaginary and 
real parts of for systems which preserve time reversal 
symmetry [40, 41, 56] : 


ELRe['k„/' 


(14) 


This parameter allows to quantify the sensitivity of 
the nth eigenmode to the openness of the system. It 
is related to the phase rigidity [pnP [37-39], where 
Pn = EL through: 


2 _ t |Pn| 

1+KI' 


(15) 


Using the above mentioned statistical assumptions, the 
spatial distribution of the normalized intensity /„ = 
|4'nP/ (|4'nP) of a complex eigenstate (where (..) cor¬ 
responds to a spatial average) reads 


R{Jn^ Pn') — 


1 

1 

1 

1_ 

Iq 

\Pn \ ^n 

-IPnP 

1 

1 

1—1 




(16) 
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This distribution interpolates between the Porter- 
Thomas distribution (|p| = 1) for closed systems and 
the Rayleigh distribution (|p| = 0) for fully open sys¬ 
tems [57]. It was formerly proposed by Pnini and 
Shapiro who considered the complex wave field as a su¬ 
perposition of standing and traveling plane waves [42]. 
This distribution had also been proposed previously by 
Zyczkowski and Lenz in the different context of the 
crossover induced by time reversal breaking in chaotic 
systems [58]. Strangely enough, to our knowledge, it was 
never shown to be valid for the complex eigenvectors of 
the non-Hermitian random Hamiltonian Since 

(hence p„) was demonstrated to be a distributed quantity 
[41, 50, 59], the distribution of the normalized intensity 
of eigenstates over the statistical ensemble is given by 

P{i) = j Pp{p)p{hp)dp. (17) 

In the weak coupling regime defined by: 

v/Var(r„)/A < 1, (18) 

and for small or moderate modal overlap d, the distri¬ 
bution of the complexness parameter of chaotic open 
wave system has been derived in Refs. [41, 59] by using 
perturbation theory, and reads 

-^(i + ^) 

Px-,M(Xn) = --- 

24Ma:2 ^1 -|- 

where Xn = 9^/d^ and M is the number of open channels. 
The distribution Pp{p) follows from Pp(p) dp = Px(x) dx 
and is given by 


2 _ f (l-|p„|) \ 

Pp-AmM\Pu\) Px-M ^ J . 

(20) 

The number of open channels M is the only param¬ 
eter of the distribution (19) and is also related to the 
fluctuations of the widths. In the weak coupling regime, 
denoting |n) the nth eigenstate of H, and by considering 
the anti-Hermitian part of Tdgff as a perturbation of H, 
the perturbation theory gives 


M 

rt, = {n\VV^\n}=J2vf. ( 21 ) 

1=1 

Thus corresponds to a sum of squares of independent 
identically distributed Gaussian random variables with 
variances cr^ [Eq. (9)]. This yields a mean width 

P" = a^M, (22) 

and a variance 

(23) 


and the rescaled widths = P^/cr^ are chi-square dis¬ 
tributed with M degrees of freedom [36, 41, 50, 59]: 


Ry;M(yn) 


.MI2-1 -'r^.12 

2^/2r(M/2) 


(24) 


Therefore the first two moments of P^ are directly related 
to M through 


M (P*^)^ 

T" Var(r(,)’ 


(25) 


and the normalized width r„ = P^/P*^, follows the dis¬ 
tribution: 

where M is the only parameter since a = b~^ = M/2. It 
is important to note that the number M and the modal 
overlap d are essential parameters for open chaotic sys¬ 
tems. Indeed, by combining the expressions (22) and 
(10), the coupling strength k can be expressed in terms 
of these parameters by: 


ttP" 

2AM 


nd 


(27) 


One can also note that Eq. (25) permits to restate the 
weak coupling condition (18) in terms of M and d as 



All the above mentioned results are checked and illus¬ 
trated in Eigs. 2 and 3. For each figure, we have diago¬ 
nalized 150 random with N = 700. We have kept 
only eigenvalues and their associated eigenvectors in the 
center of the semi-circular law (e ~ 0), where one can 
assume that A is almost constant. The results shown in 
Figs. 2 and in 3 are respectively obtained for two differ¬ 
ent values of the modal overlap and number of coupling 
channels, namely d = 0.5 and M = 20 and d = 1 and 
M = 10. The parameter d^2/M is respectively 0.16 and 
0.45 for the results associated with Figs. 2 and 3. Thus, 
we can assume a priori that the results shown in Fig. 2 
must correspond to the weak coupling regime whereas 
the results shown in Fig. 3 will deviate from this regime. 
This point is indeed supported by graphs on top where 
the empirical distributions of normalized width P„ [black 
histograms in Figs. 2(a) and 3(a)] and the empirical dis¬ 
tributions of modal phase rigidity pn [blue histograms in 
Figs. 2(b) and 3(b)] are respectively compared with the 
perturbative predictions (red continuous curves) Pr-M 
[Eq. (26)] and Pp.{M,d} [Eq. (20)], with M and d tak¬ 
ing the values used in RMT simulations (see captions). 
Only the results of Fig. 2 tally with the weak coupling 
assumption. For Figs. 3(a) and 3(b), the visible devia¬ 
tion between the histograms and the continuous curves 
is due to a larger value of the parameter dyj2/M for 


Var(r;) = 2a^M, 
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FIG. 2. (color online) The histogram corresponds to (a) the 
width distribution, (b) the modal phase rigidity distribution, 
and (c) the intensity distribution obtained by RMT calcula¬ 
tions using 

an ensemble of eigenvectors of 150 random Jfgg- with 
N = 700, M = 20 and d = 0.50. The continuous red curves 
correspond to the weak coupling predictions given by 
Eq. (26) for (a), Eq. (20) for (b) and Eq. (17) for (c) with 
Pp Pp]{M,d} • 


FIG. 3. (color online) The histogram corresponds to (a) the 
width distribution, (b) the modal phase rigidity distribution 
and (c) the intensity distribution obtained by RMT calcula¬ 
tions using an ensemble of eigenvectors of 150 random ITgfj 
with N = 700, M = 10 and d = 1. The continuous red curves 
correspond to the weak coupling predictions given by Eq. (26) 
for (a) and Eq. (20) for (b). In panel (c), the blue continuous 
curve is obtained through Eq. (17) where Pp is given by the 
discrete empirical distribution [blue histogram shown in (b)]. 


which a weak coupling assumption is clearly no longer 
valid. For completeness, we also compare the exact num¬ 
ber of coupling channels with the nearest integer given 
by expression (25) and with the value given by fitting the 
empirical distribution of with Pp;{M,d} where M is left 
as a free parameter and d fixed by RMT simulation. For 
the results corresponding to Fig. 3, there are clear dis¬ 
crepancies between the values returned by the fit M = 13 
and 2F^/Var(F„) = 7 whereas for the results correspond¬ 
ing to the weak coupling regime (Fig. 2) the estimation 
of M through expression (25) and by fitting pp are both 
consistent with M = 20. The reader’s attention should 
be drawn to the fact that the different ways of estimat¬ 
ing the number of channels [either through Eq. (25) or by 
fitting Pp] are indeed mutually consistent as long as the 
perturbative regime holds and gradually disagree when 
leaving it. In Figs. 2(c) and 3(c) we compare the empiri¬ 
cal distributions of all normalized intensities of the above 
mentioned eigenvectors (black histograms) with the for¬ 
mula (17) [red continuous curve in Fig. 2(c) and blue 
continuous curve in Fig. 3(c)]. Irrespective of the cou¬ 
pling regime, a remarkable agreement is observed over 


more than four orders of magnitude. To the best of our 
knowledge, this is the first time that the eigenvectors of 
random matrices associated to are shown to follow 
the distribution (16). However, there is a slight difference 
between the theoretical curve plotted in both figures. In¬ 
deed, the red curve of Fig. 2(c) has been obtained by sub¬ 
stituting in formula (17) Pp for the perturbative Pp;{M,d} 
[red curve shown in Fig. 2(b)], whereas because of the 
departing from the weak coupling regime, the blue curve 
of Fig. 3(c) have been built from the empirical Pp [blue 
histogram shown in Fig. 2(b)]. If in this regime, instead 
of using the empirical Pp, we use the Pp-{M,d} shown in 
Fig. 3(b), the resulting distribution differs from the blue 
curve by more than 12% near the most probable values 
and by more than 20% on the tail. In both cases, the 
Porter-Thomas distribution has the wrong tail behavior 
and deviates more than 10% from the red and blue curves 
in Figs. 2(c) and in 3(c), respectively. 
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III. VECTORIAL RESPONSE IN LOSSY EM 3D 
CAVITIES 


We now consider the EM response in a RC as a con¬ 
tinuous function of the excitation frequency /. The elec¬ 
tric field, E(r, f) created by an arbitrary current source 
Jeiro, f) localized in a volume Vq and oscillating at fre¬ 
quency f, is the solution of the vectorial Helmholtz equa¬ 
tion [60, 61]: 

V X V X E — E = —2TriffioJe, (29) 

where k = jc. The electric field can be writ¬ 
ten in term of the convolution of the source term with 
the so called Dyadic Green’s function (DGF), denoted 

Gir,fo,f), [60, 61]: 


Eir, f) = -iufio f G(f, fo, /) • Je(ro, f)dfo, (30) 
Jvo 

Indeed, the DGF is the solution of 
V X V X G(r, ro, /) - /c^G(r, tq, /) = lld(r - fo). (31) 


Since we are only interested in modeling the statistical 
properties of the EM response, the DGF appears to be an 
appropriate tool. The DGF is often written as a matrix, 


G{r,fo,f) 


Gxx 

Gyx 

G,„ 


G xy 
Gyy 
Gzy 


Gxz \ 
Gyz 

Gxz J 


(32) 


where each column (associated to the second subscript 
j = X, y, or z) stands for the three components of the 
electric field vector at the spatial coordinate r: 


( G., \ 

Ej{r,f) = G{f,foJ)-ej= Gyj . (33) 

V / 


Thus the vector {Gxj,Gyj,Gzj} is the EM response cre¬ 
ated in the cavity by a point-like elementary current 
source located at fg, polarized along the Gartesian unit 
vector f), and oscillating at frequency /. Following the 
formalism given in Refs. [31, 61], the DGF can be also 
expanded over the resonances to read: 


G(f,fo,/) 


E 


En{r) ® F„(fo) 
{kl - fc2) 


(34) 


where k = 2Trfjc. In the latter expression, losses lead 
to complex (non-real) eigenmodes (resonance states) 

[50, 56], associated with complex eigenvalues 


k 


n — 


27r/n 

c 



(35) 


the nth resonance bandwidth T^ = fn/Qn- In expression 
(34), the irrotational contribution has been omitted since 
it can be neglected far from the source region (f fg). 
The amount of overlap between adjacent resonances can 
be evaluated through the mean modal overlap, here de¬ 
fined by the ratio of the mean bandwidth T-^ and the 
mean resonant frequency spacing A-^: 

d = r^/A-''. (36) 

In the following, as we are only concerned with EM 
waves, we will omit the / superscript. 


A. Bi-orthogonality and phase rigidity for vectorial 

fields 

Equation (7) can be re-written in the spatial represen¬ 
tation as Eq. (34), which is an extension of the case of 
the scalar field, where 


G(f,f ,*) 


E 


4'n(r)4'„(f) 

U2 1.2 
fi 


which is a solution, inside of the cavity, of 
(Af + k^) G{f, 7^, k) = S{r — r*) 


(37) 


(38) 


together with the boundary conditions imposed by the 
coupling channels [56, 62]. In contrast to the Hermitian 
problem, where the residues are given by 'k„(f)'k*(f), 
here we have the residues 4'„(r’)4'„(f) due to the rela¬ 
tion (L„j = Note that different sign conventions 

are used in the quantum case [Eq. (37)] and in the elec¬ 
tromagnetic case [Eq. (34)]. We wish to draw the reader’s 
attention to the central fact that expressions (34) and (7) 
can be made formally equivalent provided that the E„’s 
are properly normalized. This normalization is achieved 
through the following transformation of the vector field: 


£n = En/\ / En ’ E^dV, 


(39) 


leading to an explicit form of the bi-orthogonal condition 
for vectorial eigenfields [formally equivalent to (5)]: 


£n ' £md\^ — 


(40) 


where Jy ...dV is the integral over the volume V of the 
cavity. Note that, in the scalar product, the standard 
conjugation is omitted due to the constraint (6). Gondi- 
tion (40) was shown to provide a smooth transition from 
the complex wave functions of an open wave system to 
the real wave functions of the corresponding closed sys¬ 
tem as 


\\£n\\^dV>l, 


where c is the velocity of light, /„ is the central frequency 
of the nth resonance and Qn is its quality factor related to 


V 


(41) 
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FIG. 4. (color online) Imaginary vs. real part of all Carte¬ 
sian components of the electric field before (a) and after (b) 
normalization by Eq. (39). 


and tends to unity for vanishing anti-Hermitian part of 
the effective Hamiltonian Tfgff. To evaluate this transi¬ 
tion, the definition of the phase rigidity for vector fields 
is extended as follows: 


jy En ■ EndV 
lyllEr^W^dv’’ 


whose modulus reads 


\Pn\ — 


/ylIfnIIW 


< I. 


(42) 


(43) 


In the following, when using the DGF given by (34) , one 
should keep in mind that En has to be replaced by £„■ 


B. Universal statistics of modes in lossy chaotic 
RCs 

As emphasized in Ref. [31], the ideal 3D chaotic cavity 
is given by a fully asymmetric room without any parallel 
walls and with defocusing parts (focusing parts can also 
be used, with restrictions regarding the centers of curva¬ 
ture [63, 64]). Therefore, using a parallelepipedic room 
to ensure a homogeneous distribution of energy is not 
an optimal option as acknowledged by the community 
of room acoustics (see Ref. [2] and references therein), 
and more recently in the EMC community [30]. Since 
we wish to address the physical situations encountered 
in reverberation chambers, we start from a parallelepi¬ 
pedic cavity and introduce simple modifications of the 
boundary. The reverberation chambers are depicted in 
detail in Fig. 1. Moving the hemisphere along the ceil¬ 
ing will provide the necessary stirring in order to gener¬ 
ate statistically independent configurations of the cavity 
(statistical ensemble). None of the caps are centered and 
both penetrate inside the cavity to a maximum length of 
15 cm. The chaotic cavity thus obtained is a 3D realiza¬ 
tion of a dispersing billiard, the well-known Sinai billiard 
[65, 66]. By such modifications, the usable volume of the 
cavity is not significantly reduced, and simultaneously. 



FIG. 5. (color online) Gomparison between empirical distri¬ 
butions of the normalized widths (black histogram) obtained 
in the cavity 1 shown in Fig. 1 where d = 0.43 and formula 
(26) with M = 2r^/Var(r„) = 11 (continuous blue curve). 
The values of y^Var^I^/A is 0.19. 


as demonstrated in Ref. [31], the cavity thus obtained 
is an excellent realisation of a chaotic RC verifying all 
expected universal statistical properties of closed chaotic 
systems. 

In the present study, we concentrate on the range of 
frequencies around 1 GHz (close to the above defined 
LUF, near the 210th eigenstate). For cavities with such 
dimensions, where losses are due only to the finite con¬ 
ductivity of the walls, the mean quality factor Q would 
be of the order of 10^ — 10®. However, in practice, the 
latter is rather of the order of 10® due to leakage through 
antennas or defects in the walls or sundry objects un¬ 
der test. For the sake of simplicity, in our model we 
mimic the reduced quality factor by introducing sub¬ 
surfaces (colored patches in Fig. 1) of reduced conduc¬ 
tivity distributed on the walls. The total area of 

these sub-surfaces is The conductivity of copper 

{(TQn = 5.7 X 10^ S/m) is attributed to the remaining 
surface (in gray in Fig. 1). The values of are chosen 
so that Q takes on realistic values comprising between 
1500 and 2000. While modifying (passing from the 
configuration with the three red patches, now indicated 
as cavity 1, to the configuration with six patches, indi¬ 
cated as cavity 2), Q is approximately kept constant by 
adjusting ex [7, H, 67]. In our model, typ¬ 

ical values of o-^j^g are ~ 10“'^cr/.„. With such a huge 
ratio, only the non-copper parts of the walls contribute 
significantly to the losses. The eigenvalues and eigen- 
modes were obtained by solving the vectorial Helmholtz 
equation with impedance boundary conditions with the 
help of the commercial software Comsol. We noticed that 
the solver gives an arbitrary global phase to the vecto¬ 
rial electric field En associated with the nth resonance 
























[Fig. 4(a)]. The phase is identical for the three Carte¬ 
sian components of the electric field but changes from 
mode to mode. The effect of the normalization (39) is 
shown in Fig. 4(b), where we can check that the real and 
imaginary parts of Cartesian components of viewed 
as spatially distributed random variables, are thus made 
independent. The lossy patches give rise to the W term 
in the effective Hamiltonian as they are opening the sys¬ 
tem to the ‘environment’ (for details see Refs. [56, 62]). 

For 30 different configurations obtained by moving the 
hemisphere along the ceiling, we computed 50 eigen- 
modes around the 200th resonance . For a modal overlap 
d ~ 0.45, the values of Y^Var(r„)/A are of the order of 
0.15 for both cavities, so that a weak coupling assumption 
is reasonable. Indeed, a good agreement between the em¬ 
pirical distributions of the widths and formula (26) with 
M = 2r^/Var(F„) is obtained (see Fig. 5). In Fig. 5 
only results concerning cavity 1 are shown, those con¬ 
cerning cavity 2 lead to the same conclusion. Indeed 
in this cavity, the parameter Y^Var(r„)/A = 0.13 and 
the distribution of the widths is in perfect agreement 
with the distribution (26) where M is evaluated through 
M = 2rVVar(r„) = 23. 

The next step is to recover the distribution of the phase 
rigidity. It will now be shown that the relation (15) be¬ 
tween pn and permits to properly define an extended 
complexness parameter which, for the case of vectorial 
fields in chaotic cavities, where the modes are expected 
to be homogenously distributed and isotropic, must fol¬ 
low the same distribution as in the scalar case. Indeed, 
by combining Eqs. (15) and (43), one has 



The bi-orthogonality conditions (40) yield 
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FIG. 6. (color online) In cavity 1, (a) the distribution of the 
phase rigidity (blue histogram) of the 50 x 30 eigenfields is 
compared to the analytical perturbative Pp;{M,d} (red curve) 
with d=0.43 and M = 11. (b) The distribution of the nor¬ 
malized intensity (black histogram) of the three components 
of the same set of eigenfields is compared to expression (17) 
(red curve) using the analytical Pp-{M,d} shown in panel (a). 


and one can deduce from Eq. (44) the following extended 
definition for the complexness parameter : 


^Im 

f 

■ Im 

f 

) 

(Re 

f 

• Re 

f 

*-'71 

) 


(45) 


Since we are interested in chaotic cavities, we can assume 
that the eigenfields £„ are statistically isotropic so that 
the Cartesian components £n,x,£n,y,£n,z are independent 
and identically distributed random variables yielding 


2 _ (Im[£„,„]’^+Im[£„,„]’^-i-Im[£„,,]’^) 
“ (Re[£:„,,,]’'-KRe[£:„,„]=-sRe[£„,,]=) 
(Imfc/]^) 

= (RhifTp) ’ 


(46) 


where U follows the same distribution as the scalar eigen¬ 
fields rkn. Hence the extended complexness parameter 
(and the corresponding phase rigidity), as well as the 
normalized intensities of the Cartesian components of the 
eigenfields, must obey the same distribution laws as for 
the scalar case [Eqs. (19) and (17)]. This is illustrated 
in Fig. 6 for cavity 1, where the distribution of Pn (blue 
histogram) is computed for all the eigenfields mentioned 
above and compared to the perturbative Pp-{M,d} (red 
curve) with d = 0.43 and M = 11 as deduced from the 
nearest integer of formula (25). In the bottom part of 
the figure, the distribution of the normalized intensity 
(histogram) of the three components of all eigenfields is 
compared to expression (17) (red curve) with the above 
mentioned analytical Pp;{M,d}- In this case where the 
weak coupling limit is ensured, the agreement between 
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dom variables whose ratio of their mean values is fixed by 
the value of the complexness parameter [see Eq. (45)]. 
Consequently, for a given eigenfield the distribution 
of its normalized intensities is given by [68] 


Pi{In]p) 


9 In 

— , =exp 

V\p\^-\p\^ 


3J„ \ / 3I„|p| \ 

1-HV 'U-HV 

(49) 


with Ii the modified Bessel function of the first kind and 
p = Pn the phase rigidity of Note that the difference 
between this equation and Eq. (16) stems from the dif¬ 
ferent degrees of freedom in the statistical components. 
Therefore, the distribution of the normalized modulus 
squared of an ensemble of eigenfields follows 


P\\E\\^ix)= J Pp{p)Pxi^{x);p)dp, (50) 


FIG. 7. (color online) The distribution of the normalized 
modulus squared of the 1500 computed eigenfields of cavity 1 
(black histogram) is compared to expression (50) (red curve) 
with the analytical shown in Fig. 6(a). 


the empirical distributions and the RMT predictions is 
remarkable. 

We have thus demonstrated that, thanks to an appro¬ 
priate definition of the phase rigidity, the Cartesian com¬ 
ponents of eigenfields of vectorial chaotic systems and 
eigenvectors of random are statistically equivalent. 
However, some typical quantities related with a vectorial 
field problem such as the distribution of the normalized 
modulus squared of the vector field, defined for the eigen¬ 
field £n by 


T — 
-^n — 


C C* 

^n,x^n,x 


^ri,y£n,y 


c c* 


F ■ F* 

*~'r. 




^n,y^n,y 


■ F F* 
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(47) 

have no direct equivalent in RMT formalism. Neverthe¬ 
less, the knowledge of statistical properties of Cartesian 
components brings the key elements to derive the dis¬ 
tribution of In- Indeed, the numerator of (47) can be 
written as 


Fn ■ FI = Re 


• Re 


Im 


• Im 


(48) 


and, as shown above, in chaotic RCs due to the statistical 


uniformity and isotropy of the EM field. Re 
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Un 


•Re 
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Un 


and Im 



•Im 
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are two independent random vari¬ 


ables, both are equivalent to the sum of the square 
of three identically distributed Gaussian random vari¬ 
ables. But the Gaussian random variables correspond¬ 


ing 


to Re 


Re 


and those corresponding to 


Im 


■ Im 

f 

are not necessarily the same. Hence 

Re 

f 

• Re 

f 

and Im 

f 

■ Im 

f 

are two Xs 


where Pp = Pp-^{M;d) in the case of the weak coupling 
regime. The reader’s attention should be drawn to the 
fact when jpj tends to I, Px —>■ xi (which is expected 
for lossless chaotic RCs) and when p tend to 0, Px —>■ 
X§ (which is expected for fully open chaotic RCs). The 
validity of the above derivation is verified in Fig. 7 where 
the distribution of of the previously mentioned set of 
eigenfields of cavity 1 (black histogram) is compared to 
distribution (50) (red curve), where Pp is replaced by the 
perturbative Pp;{M;d} shown in Fig. 6(a). 


C. Effectiveness of random matrix to predict 
statistics of the response in chaotic RCs 

In this last section, we are interested in the distribution 
of the normalized intensity of Cartesian components or 
equivalently in the distribution of the normalized mod¬ 
ulus squared of an ensemble of EM responses in chaotic 
RCs. Originally, the distribution (16) was proposed by 
Pnini and Shapiro to describe the statistics of the wave 
function of open disordered or chaotic systems [42]. For 
a given configuration and at a given excitation frequency, 
the phase rigidity of the wave function is thus required. 
In the case of well isolated resonances, the phase rigidity 
of the wave function given by the response (7) is deter¬ 
mined by the values of the phase rigidity of the individ¬ 
ual resonance state and is really meaningful at the cen¬ 
tral frequencies of the resonances. Obviously, when the 
modal overlap comes into play, the phase rigidity of the 
response is built upon several different states which con¬ 
tribute according to their spatial overlaps with the source 
and to how far their central frequencies are from the exci¬ 
tation frequency. Therefore the distribution of the phase 
rigidity of the response is generally different from the cor¬ 
responding quantity for the individual eigenfunctions, as 
already emphasized in Ref. [39] . Consequently, for an en¬ 
semble of EM responses {E}, the distribution of the nor¬ 
malized intensity of each Cartesian component and the 
distribution of the normalized intensity of the response 
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FIG. 8. (color online) Comparison of the distribution of the 
normalized modulus squared of an ensemble of EM responses 
of cavity 2 (black histogram) and the distribution (52) where 
Pp is the synthetic distribution obtained via RMT response 
simulation with M = 23 and d = 0.43 (green continuous 
curve). For the sake of comparison, the xl distribution, ex¬ 
pected for closed systems, and Xg, expected for fully open sys¬ 
tems, are represented by the dashed and dotted violet lines, 
respectively. 

are expected to be given by 

P{I) = J Pp{p)P{I;p)dp, (51) 

P|l^l|.(x) = I Pp{p)Px{I{x);p)dp, (52) 

respectively, where p = (^E ■ p'j / ^ and Pp is, to 

our knowledge, not known analytically and only empiri¬ 
cally accessible through a statistical sample of values of 
p. However, if the number of effective coupling chan¬ 
nels M and their coupling strength k are known, the 
statistical behavior of any open chaotic system can be 
simulated through random Pgff simulations. Especially, 
the distribution of the phase rigidity of an ensemble of 
responses of a realistic chaotic system should correspond 
with the phase rigidity distribution of an ensemble of 
random matrix responses (7) built upon the eigenvec¬ 
tors and eigenvalues of random Pgff, where M and k are 
those of the corresponding realistic system. This syn¬ 
thetic distribution of p produced by RMT simulations 
can then replace Pp in Eqs. (51) and (52) and the result¬ 
ing distributions should predict respectively the intensity 
distribution of Cartesian components and the distribu¬ 
tion of normalized modulus square of an ensemble of EM 
responses of chaotic RCs. We now demonstrate the rele¬ 
vance of the above surmise in the case of an ensemble of 
response computed in the chaotic electromagnetic rever¬ 
beration chamber cavity 2. This ensemble is obtained by 


varying the configuration given here by the threesome: 
excitation frequency, polarization of the source and po¬ 
sition of the hemisphere. More precisely, we computed 
50 eigenmodes of cavity 2 around the 200th resonance 
for 30 positions of the hemisphere, and built the re¬ 
sponse using Eq. (33) (evaluated at 16 different points 
inside the cavity) by means of the DGF (34), for 300 
equidistant excitation frequencies comprised in the inter¬ 
val [969.4 MHz, 1024 MHz]. For the results presented be¬ 
low, the average modal overlap is d = 0.43 corresponding 
to a mean quality factor Q = 1500. As mentioned above, 
in this frequency range, a weak coupling regime can be 
assumed for cavity 2. Thus the number of coupling chan¬ 
nels as well as the coupling strength can be estimated 
thanks to Eq. (25) M = 2r^/Var(r„) = 23 and Eq. (27) 
K = T:d/(2M) = 0.03. In Fig. 8, the distribution of the 
normalized modulus square of the above mentioned en¬ 
semble of EM responses is compared with the distribution 
(52) where Pp is the synthetic distribution obtained via 
RMT response simulation with M = 23 and d = 0.43. A 
good agreement between histogram and RMT prediction 
is observed. In contrast, deviations from the predictions 
for the closed and fully open systems are found. Thus we 
find universal statistics for all fully-chaotic partly-open 
reverberation chambers with the same number of effec¬ 
tive channels and modal overlap. 


IV. CONCLUSION 

Starting from the effective Hamiltonian i7gff formal¬ 
ism including random matrix theory, we can derive the 
distribution of intensities for scalar waves in the weak 
coupling limit. The only parameters entering our model 
are the modal overlap and the number of open chan¬ 
nels. For weak coupling (corresponding to small or mod¬ 
erate modal overlap) we hnd good agreement between 
random matrix eigenvectors and our theoretical predic¬ 
tions, whereas, in the range of strong coupling, deviations 
are found. The discrepancies are essentially due to the 
deviations in the distribution of the phase rigidity. For 
further comparison, we calculated the electric field dis¬ 
tribution inside a lossy chaotic reverberation chamber. 
By comparing the thus obtained width distributions with 
the x^ distribution expected in the weak coupling regime, 
we extract the only two parameters, namely the modal 
overlap and the number of open channels. Concerning 
the vectorial electric eigenfield, a good agreement is seen 
for the distributions of (i) the phase rigidity |p„|, (ii) the 
individual Cartesian component intensities /„ as well as 
(iii) the held intensity I. Finally, we derive the distribu¬ 
tion for the response, typically the quantity of interest 
in reverberation chambers used for electromagnetic com¬ 
patibility. The obtained distribution corresponds well 
with the theoretical one, where the synthetic distribu¬ 
tion of the phase rigidity obtained from random matrix 
numerics is used. We would like to emphasize that, in 
all cases, the only parameters are the modal overlap and 
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the number of channels, which are known a priori. 
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